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Abstract 

We study the final state strong interactions of the B DK decay pro- 
cesses, using the Regge model. We conclude that the final state interaction 
phases are very small, typically a few degrees. Neglecting final state interac- 
tions in obtaining the weak decay amplitudes is a good approximation. 
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The importance of studying the final state strong interactions of the two body 
non-leptonic B decays is due to the desire of determining the relevant CKM matrix 
elements and of studying the CP violation effects in the B decay processes. In the 
first case the final state strong interactions will modify the value of the CKM matrix 
elements extracted from the bare weak decay amplitudes. In the latter case, in order 
to have observable CP violation effects, it is necessary to find the difference between 
the decay rate of a process i ^ f and the rate of its charge conjugate process, 
i f where there have to be at least two interfering amplitudes with different weak 
interaction phases as well as different strong interactions phases. By writing the 
bare weak interaction T matrix as, 

A = g^Ai + ^2^2 , (1) 

where the Qi are the weak couplings, we have, after taking the final state strong 
interactions into account, 

</|A|z>=(7iAie*-^+^72v42e^"^ , (2) 

< f\A\i >= glA.e''^- + glA^e''^- , (3) 

where the are the strong interaction phases. The difference in the decay widths 
of i — s> / and i ^ f is, 

r - f ~ Im{glg2) sin(ai - ^2) . (4) 

We see that whether there are observable CP violating effects depends crucially on 
the strong interaction phase difference between the two interfering amplitudes. 

Despite of the importance of the final state strong interactions in B decay pro- 
cesses, it is however difficult to deal with and poorly understood. In the K decay 
system the final state strong interaction phase difference can be estimated reliably 
using the low energy effective theory of strong interactions. For the D decay system 
in which the center of mass energy is too large to apply the low energy effective 
theory, a single resonance dominance assumption was used in studying the final 
state strong interactions. Under this assumption, the strong interaction amplitude 
is simply parametrized as the function of the decay rate of the resonance and the 
couplings between the resonance and the interacting particles. The phase difference 
can therefore be estimated in a resonable range of parameter space [|| . This method 
can not be successful for the B decays because we know from QCD that the number 
density of resonances is an increasing function of the mass. 

In this paper, we study the final state strong interactions of the two body non- 
leptonic B decay processes using the Regge model analysis 0. This method is based 
brief discussion based on the Regge model was given in ref. B 
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upon the strong interaction duality argument: When the center of mass energy gets 
high enough, the summation of the contributions from the s-channel resonances to 
the amphtude is equivalent to the summation of the contributions from leading t- 
channel Reggeon exchanges. The Regge model has been proven to be very successful 
in explaining high energy hadron scattering in the small — |t| region 0. As an 
example, in the case of vrA^ elastic scattering resonance models only work when 
the center of mass energy is less than s ~ 6GeV^ and fail badly when s exceeds 
lOGeV"^ where the Regge model starts to work very well [|] (In charge exchange 
processes it works even at lower energies). In order to have a clear insight on the 
physics concerned, we limit ourselves to the processes without penguin diagram 
contributions. The absorbtion part of the amplitudes in such a case come only 
from "soft" final state strong interactions. Especially, we discuss the B — > DK and 
B DK processes although applications of our method to any two body decay 
process are straightforward. 

The relation between the full physical amplitude when the final state interactions 
are taken into account and the bare amplitude is given by the final state theorem 
of Migdal and Watson 0, 

Atot = S^Abare ■ (5) 

where S is the final state strong interaction S matrix for the given partial wave 
(J = in the present case). 

The final states are classified according to their flavor quantum numbers. For 
a fixed final state f various intermediate states can contribute coherently to the 
rescattering. For example for the final state D^K~ both D^K~ and D^K^ are 
possible as intermediate states leading to a 2 x 2 S* matrix. While for D^K' , D~K^ 
case things are more complicated since the D~ri and D~7i^ also contribute and the 
S matrix is 4 x 4 ^. The Pomeron and Reggeon contributions to the s channel 
scattering amplitude A(s,t) can be written down systematically using the crossing 
symmetry between the s channel and the t channel amplitudes, the line-reversal 



law P|, and the SU{3) relations P], 

AiD'^K- ^ D'^K') = P + p + A2 + u; + f, (6) 

A{D^K^ _ D^K^) = P _ p _ ^2 + ^ + / , (7) 

A{D-K^ ^ D-R'') = P + p + A2 + u; + f, (8) 

AiD+K'^ D+K^) =P-p + A2-uj + f, (9) 

A{D^K- ^ D^K-) = P + p-A2-uj + f, (10) 



■^There are many other states with the same total spin and flavor quantum numbers, includ- 
ing two body states with excited particles and multi-particle states. As an approximation we 
systematically neglect them. We will discuss their influence later in this paper. 
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(11) 

(12) 
(13) 
(14) 
(15) 

where the P, p, u;, A2, f and /' denote the contribution from the Pomeron, p, u, 
A2, f and /' Reggeon exchanges to the the amphtude, respectively. For the charge 
exchange processes, we have 
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where the K* and K** denote the K* and K** Reggeon contributions, respectively. 
The SU{3) coefficients are written exphcitly in above formulas. Please notice that 
the Pomeron couphngs may depend on the different flavour content while various 
Reggeon contributions are related to each other by the strong exchange degeneracy 
(SED) requirement 0. Especially, the magnitude of each Reggeon contribution in 
an amplitude should be equal. The SED is known to hold quite well for vector and 
tensor Reggeon exchange processes especially when the strange quark is involved. 
Sign differences between Reggeon contributions in the above formulas lead to the 
cancellation between the imaginary part of the Reggeon contributions to s channel 
exotic amplitudes. This is just a manifestation of the fact that the absence of an 
imaginary part in A^ implies that there is no s channel resonance, because of the 
duality argument. 

We parametrize the Pomeron and the Reggeon exchange amplitudes in the small 
— \t\ region in the following way: 

/3^(t) ( - ) e-^i°"W , (23) 



P 

' '.So. 



R = (-) . n ' (24) 



(3P{t) = pPe''''' and /3^(t) = /3^e'^'^* , (25) 
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where ± sign refers to odd/even signatures of the exchanging Reggeons. The SU{3) 
invariant couphngs (3^ are normahzed such that it is just the value of the residue 
function of the Reggeon-matter couphng at origin, i.e., t = 0. In accordance with 
the traditional treatment we use an exponential form of parametrization to charac- 
terize the form factor of the Reggeon-matter couplings. The parameter can be 
estimated using the Veneziano model for meson-meson scatterings. 



which obeys approximately the relation: = where mji is the mass of the 
corresponding elementary particle of the Reggeon R. We will approximate by 
l/m|. in the following. 

It is worth pointing out that our parametrization, eq. ([2^), is different from what 
is usually used in the literature: 

A\s, t) = -a;/3^e" * (^yj g"*™ ('^/^ (27) 

with X = 1/ — i for even/odd signatures, respectively. The reason is that as 7^ 0.5, 
the parametrization eq.(p7|) violates SED by introducing an imaginary part to the 
Regge amplitude for a s-channel exotic process and is therefore not adequate for 
the present discussion. For the trajectory functions a^(t) and a^(t) we use the 
following values P, ^, 

a^(t) ^ 1.08 + 0.25t , (28) 

aP{t) c^0.U + 0.93t , (29) 

a^*(t) ~ 0.3 + 0.93t , (30) 

a/'(t) = 0.1 + 0.93t . (31) 

To estimate the residue functions l3^{t) and /5^(t) we assume factorization. That 
is, for a process AB CD we have, 

(3AB-,CD{t) = (3Ac{t)f3BD{t) , (32) 

both for f3^{t) and I3^{t). In the literature, the Pomeron coupling to the pion is 
taken as a single pole form in the small —\t\ region |§, 

The t dependence of Pomeron couplings to D and K mesons are not known experi- 
mentally, we simply assume they are the same as that in eq.(^)0. For the residue 

"^Roughly speaking, the formfactor is the Fourier transformed version of the size of the hadron. 
Since the tt, D and K mesons both involve hght quarks their size is expected to be similar. 
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function of the Pomeron amplitude f3^{t), we take 

r (t) = « - - { i-t/o.7i ) ' " = • (^^^ 

The value of P^j^ can be estimated as follows. From up and /^p data one gets, 

pP^su) ^ ^P^'iuu) , (35) 

and 

/3^(cw) ^ , (36) 

from [^. Using isospin invariance 

/3^(mm) = (3^{ud) = (3^{dd) , (37) 

and the additive quark counting rule we obtain 

f3^{DK) = (3^{cu) + (3^{cs) + (3^{us) + (3^{uu) ~ 12.0 . (38) 

The coupling constant in eq. (^) cannot be estimated from SU{3) constraints, 
since the D meson is a SU{3) triplet and may have different coupling comparing 
with the SU{3) octet mesons. It can however be fixed from the SU (4) relation 

f3^ = (3§, (39) 

provided that the coupling constant can be estimated from A^A^, vrA^ and KN 
scattering data. Eq.(^) is also a direct consequence of the vector meson coupling 
universality, which is respected satisfactorily in the 7t{K)N processes. 

The Regge pole amplitudes for 'k{K)N elastic scattering can be written down 
analogously. We have fitted the experimental data on the t[{K)N total cross-sections 
and find that ~ 3.11. 

At this stage we are able to calculate the two body scattering Regge amplitudes 
given above which is however not directly applicable for studying the final state 
interactions of the B decay system. Because the B meson has spin 0, the corre- 
sponding scattering amplitude should be the s-partial wave projection of the full 
Regge amplitude: 

^j=o(s) = A{s,t)dt, (40) 

where 

A = A(s, m^, m\) = + + — 2sm\ — 2sml — 2m\m\ , (41) 

and 
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Before performing the numerical calculation we first remark on the reliability 
of our method. It is known that the Regge parameterization is only valid in the 
small —\t\ region. In order to obtain the s-wave amplitude we need to perform 
the integration over the whole - \t\ region, and further, low partial wave projections 
of the Regge amplitudes are considered as the most unreliable part of the theory. 
However, the large — \t\ contributions to the total scattering amplitude is suppressed 
by 1/s, therefore the uncertainty due to the invalidity of the Regge amplitude at 
large — \t\ region and the dependence of our results on the parameterization form in 
the small —\t\ region should not be important. For the large —\t\ region Regge cuts 
and exchange channel Reggeons dominate. The Regge trajectory in the u-channel 
are D* and D**(2460) in our case. However, these u-channel Reggeons have a very 
small intercept a^j: 

2 

aO,. = al„ = 1 - , ""^^ , ~ -1 , (43) 

their contributions are therefore negligible. 

For the Regge cuts, it was realized long ago that one of the main defects of Regge 
poles is that they give too large a contribution to the low partial waves. Some extra 
absorption provided by the Regge cuts is necessary. However, the absorption effects 
in our case (meson-meson scattering) should not be large either. There is strong 
experimental evidence indicating that the absorption effects are much smaller in 
T^N scattering than in A^A^ scattering. It is natural to expect that in meson meson 
scattering these absorption effects are even smaller. This may be understood from 
the simple Reggeized absorption model |]T0[ despite of the fact that Regge cuts are 
less well known theoretically than Regge poles. In the simple Reggeized absorption 
model the cuts generated by two Pomerons or one Pomeron one Reggeon exchange 
contribute to the full s-partial wave amplitude as: 

^p+p{g)p = ^p[^i__M!_fAy] , (44) 

and 

where e = 0.08, cp = + a'p (\og{-^) — TT/2i^. The phenomenological enhancement 
factor A characterizes the contribution from quasi-elastic intermediate states. We 
read off from eqs.(|^), (|45|) that the contribution from the cuts is proportional to 
which is much smaller in our case, ~ 12, than in the case of vrA^ scattering. 
This is crucial to reduce significantly the correction to the s wave amplitude, different 
from the ttN scattering case. For example, taking A/j = 2 in the latter case will lead 
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to a complete absorption in the s-wave amplitude, As ^ while the same value 
of the Xji parameter only reduces As{DK) by ~ 30%. Even though we have poor 
knowledge in estimating the Xr parameter we expect this amount of reduction is 
resonable, as an upper limit. The Xp parameter may be safely neglected. According 
to ref. [|] the inclusion of the P{g)P cut in the A^A^ scattering case only leads to a 
small change of the f3 parameter. 

Because of the above arguments we use in the following only the simple Regge pole 
model to evaluate the effects of the final state interactions of the B decay processes 
from the final state theorem, eq.@. In the practical calculation it is convenient to 
study the problem in the strong interaction eigenstates for which the 5* matrix has 
a diagonal form. For the DK case we have 

where Ai^o = P + 2{p- A2), Aj^i = P + 2(^2 - p) and \DK >j=o= j={\D+K- > 
>, \DK >/=!= -j^ilD+K- > +\D^K^ >. For the DK case, there are two 
degenerate states of each isospin which are, 

|0^ >= — ^{\D°K~ > +\D-K° >), 10^ >= \D;r] >, (47) 
v2 

|1^ >= — ^{\D°K- > -\D-K^ >), 11^ >= |D;7r° > . (48) 
V2 

The 7 = 05" matrices, Sij = Sij + piPjAf~^ {pi = ^/Xi/s) are the following: 

(49) 

In the strong interaction eigenstates 5* is parametrized as, 

S^=° = diag{r]ie^'^'), (50) 

The Tji parameter characterizes the inelasticity of the given process. 
The numerical results for the parameters given above are|, 

Vo{DK) = 0.81, So = -3.3°; rii{DK) = 0.84, Si = 1.8°; (51) 



10^ > 
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|02>\ 
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<l'l 


ll'>\ 


\0' > 


<02 


|02 > j ' 




111 > 


< 1^1 


ll'>/ 



^The S matrix is equal to, 

' a b 



b a' 

with a — a' ,b ^ O(10^^)a. The S^^^ matrix is of the foUowing form: 

X y 
y x' 

with x,x' >> y. Approximately x = ^/a, y — b/2^/a and x' = \/a+ (a' — a)/2^/a. The numerical 
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Vo{D;r]) = 0.85, So 



lT;r]i{D;n^) = 0.83, (5i 



-1.5°; 



(52) 



r}o{DK) = 0.83, So = 2.2°;r}i{DK) = 0.83 ,Si 



-1.8°. 



(53) 



Prom above we see that the phase shifts S are very small (modulo utt, of course). 
This can be understood qualitatively: At the center of mass energy niB the Pomeron 
contribution is dominant which gives almost a purely imaginary contribution to the 
A amplitude and therefore, the S matrix elements are almost pm^ely real. Purther- 
more, in our case the Pomeron contribution is much smaller than in the vrA^ or NN 
case because of the smaller couplings. Only small cancellation occur in the S matrix 
elements between the Pomeron contribution (mainly a negative real value) and unity 
(the value of the S matrix element in the limit of vanishing interaction). 

We notice that in our results the rj parameters are slightly less than 1, which 
indicates that we have neglected some final states in addition to those considered. 
In order to restore unitarity one has to take them into account which will lead to 
corrections to the relations between the bare amplitudes and the full amplitudes 
obtained above. The calculation of the effects brought by these states go beyond 
the ability of the present model analysis. In the following we try to give a simple 
qualitative estimate. 

Consider a given strong interaction eigenstate \a > which diagonalizes the S 
matrix in the incomplete basis like the states that we have considered: < a\S\a >— 
776^**1 and 77 < 1. Since rj is less than 1 unitarity tells us that in the complete 
basis the S matrix is still nondiagonal. Many new states can contribute to restore 
unitarity by introducing the non-diagonal matrix element < new\S\a >. Now we 
assume that all the effects of these states are characterized by a single effective state 
I a' > which, together with \a > leads to a 2 x 2 unitary S matrix which can be 
simply parametrized as. 




(54) 



The solution of S^/^ is 




(55) 



results are, 





0.91 + 3.06i- 10-3 
2.76i • 10-2 



2.76z • 10-2 
0.93 + 2.52M0-3 



We see that the nondiagonal elements are much smaller than the diagonal one. 
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where 

I + 2r] cos(5i — ^2) + v"^ 



X 



\ 2{l + T]cosi6i-62)) 



(56) 



and 



^1 + ^2 , . -1 r/sin((5i - , , . 

ai2 = ± sm . (57) 

The full physical amplitude we are interested in < a\B >^ can then be written as, 

< a\B >f=xe''^ < a\B +Vl - xHe'^"-^^ < a'\B >'' . (58) 

From eqs. (|56|) and (^) we see that x is very close to 1 and ai very close to 61 when 
7] ~ 0.8 within a reasonable range of the 62 parameter, say, \62\ < 30° (remember 
that 62 should also not be very large because of the Pomeron dominance). Therefore 
the phase shift ai remains small. This is in agreement with the present experimental 
evidence. Furthermore since < a'\B > represents an averaged value of contributions 
from many states it is expected that large cancelation between different amplitudes 
should occur which leads to| < a'\B > / < a\B > \ <<1. Therefore we can 
conclude from eq. (|SHD that the bare weak decay amplitude is a good approximation 
to the full physical amplitude, accurate up to, roughly speaking, about 10 percent. 

To conclude, the partial wave elastic unitarity is a good approximation for the 
DK scatterings at s = even though we know that " Pomeron dominance" implies 
inelasticity. In our special case it happens that these two statements are consistent. 
One simple way to understand qualitatively the smallness of A6j (= 61 — 6q) is that 
only the Reggeon contribution is isospin dependent. Therefore it is proportional 
to Ai?//(1 — P) and be a small quantity. The above analysis can be applied to 
other processes. In Dn processes similar qualitative results should hold. In the vrvr 
and Kit cases the Pomeron contribution gets larger (by a factor of ~ 2), the value 
of the rj parameter is further decreased (r/ ~ 0.6 — 0.7). Uncertainties raised from 
inelasticity are more serious and sound conclusions are difficult to make. 
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